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ABSTRACT 

The effect of a sharp front separating the quark-gluon plasma phase from the 
hadronic phase is investigated. Energy- momentum conservation and baryon number 
conservation constrain the possible temperature jump across the front. If one assumes 
that the temperature in the hadronic phase is T ^ 200 MeV , as has been suggested by 
numerous results from relativistic ion collisions, one can determine the corresponding 
temperature in the quark phase with the help of continuity equations across the front. 
The calculations reveal that the quark phase must be in a strongly supercooled state. 
The stability of this solution with respect to minor modifications is investigated. In 
particular the effect of an admixture of hadronic matter in the quark phase (e.g. in the 
form of bubbles) is considered in detail. In the absence of admixture the transition pro- 
ceeds via a detonation transition and is accompanied by a substantial super-cooling of 
the quark-gluon plasma phase. The detonation is accompanied by less supercooling if a 
small fraction of bubbles is allowed. By increasing the fraction of bubbles the supercool- 
ing becomes weaker and eventually the transition proceeds via a smoother deflagration 
wave. 
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1. Introduction 



Experimental results from heavy-ion collisions at CERN indicate that hadrons are 
being produced with a temperature of about T ~ 200 MeV [1,2,3,4]. From an analysis 
of the ratios of strange baryons as measured by the WA85 collaboration at CERN 
[2] one finds that a baryonic chemical potential oi jiB = 300 MeV is consistent with 
their results [5-12] (provided no transverse flow exists [13]). As a first attempt one can 
therefore assume that the thermodynamic properties of the hadronic gas are known. It 
is then of interest to investigate whether this hadronic gas could have been in a quark- 
gluon plasma phase before freeze-out. To this end we assume that a sharp front existed 
separating the two phases. Various conservation laws , e.g. baryon number, energy 
and momentum, restrict the temperature difference across the front. In addition, the 
second law of thermodynamics requires the entropy to increase across the front. It 
can easily be established that continuity in the various conservation laws leads to only 
one possibility, namely a highly supercooled quark-gluon plasma must have existed 
before the transition to the hadronic phase. The degree of supercooling is surprisingly 
high, Tq ~ 70 MeV. This leads one to suspect whether some alternate mechanism 
started operating, preventing the quark-gluon plasma from supercooling to such a large 
degree. In this paper we would like to investigate several possibilities which could 
prevent this. In particular, it is known that supercooling enhances the formation of 
bubbles of hadronic matter. We will therefore consider that the quark phase was not 
pure but that some admixture of hadronic matter exists. It is known that the boiling 
of liquids only takes place if impurities are present or if some admixture of the gaseous 
state exists. Our results confirm this standard knowledge. The transition from quark- 
gluon plasma to hadrons proceeds in a violent manner towards superheated hadrons. 
The presence of an admixture of hadronic material in the form of bubbles or in any 
other form makes the transition easier. Our results agree with this expectation, from 
the moment the fraction of hadronic material is appreciable, the violent detonation 
type transition turns into a smoother deflagration type. The transition between the two 
solution shows an S-type of curve. This means that for a given reaction the solution is 
not unique. A possible interpretation would be that in this region a chaotic turbulent 
type of transition takes place. 

The paper is organized as follows : in the next section we discuss the equations 
of state we use in the the quark-gluon plasma and in the hadronic phases; in section 
3 we present for completeness the equations expressing energy, momentum and baryon 
number conservation. We also discuss the adiabaticity condition; in section 4 we discuss 
the admixture of "impurities" into the quark phase in the form of hadronic matter. The 
velocities around the sharp front are discussed in section 5. In section 6 we present our 
conclusions and give a summary of our results. 

2. Equation of State in the Quark Gluon Plasma 

In this section we discuss the equation of state used in the quark- gluon plasma 
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phase. We consider the cut-ofF model, motivated by numerical results from lattice gauge 
theories [14]. The precise value of the cut-ofF has a big influence on our results. For 
comparison we also consider the M.I.T. bag equation of state which has been previously 
discussed in detail [15]. In the cut-ofF model one leaves out the low momentum part of 
the particles, e.g. for gluons we have : 



It is well-known that the cut-off model reproduces the lattice results reasonably 
well. It is not excluded that a better way to proceed is to give the gluon a temperature 
dependent mass as has been argued recently by Goloviznin and Satz [16]. This proposal 
cannot be followed up here because we want to include the possibility of super-cooling 
and heating and we, therefore, have to know the equation of state both below and above 
the critical temperature. We have chosen a value of 5Tc for the cut-off parameter in 
order to reproduce the critical temperature Tc obtained in lattice gauge theories for 
A*B = [17]. Certain quantities are very sensitive to the precise value of this cut-ofF and 
to show this influence we also present for comparison results obtained with a value of 
3.5Tc. The corresponding results for the energy density are shown in Fig. 1. On the 
hadronic side we incorporate the hard core radius of hadrons by choosing [18] 



where the subscript zero refers to point-like quantities, Vq refers to a typical proper 
volume of a hadron which we take to have a radius given by 0.8 fm. The hadronic side 
in our approach is a composition of non-interacting quantum gases. The use of quantum 
statistics is essential since in a wide range of fis and T the Boltzmann approximation 
is not satisfied. We include all well established hadronic resonances listed in the latest 
edition of the Review of Particle Properties [19]. A correction similar to equation (11) is 
performed to all thermodynamic quantities. The behaviour suggested by the M.I.T.-bag 
equation of state is very different and is not supported by the lattice results. The larger 
cut-off reduces the energy-density in the quark-gluon plasma since a big contribution is 
being left out. The latent heat thus becomes very small for a large cut-off. Since lattice 
results are only available for zero chemical potential, fiB=0, it is not immediately clear 
how to generalize the cut-off model to finite baryon densities. Assuming the cut-off to 
be independent of hb we have determined the phase diagram depicted in Fig. 2 where 
we compare it with the M.I.T.-bag model one. At //s = the two are chosen to coincide. 
For non-zero values of jib they become more and more different. For the region we are 
interested in //s ~ 300 MeV and the difference is not very large. Having discussed the 
equations of state we will be using in each phase, we now proceed to the equations of 
continuity across the sharp front. 




(1) 
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3. Equations of Continuity across the Boundary 



Energy and momentum conservation across the front leads to the weU-known two 
equations (see e.g. [20]) : 

{eh + Ph)llvh + Ph= {eq + PqhqVq + Pq, (3) 

(en + Ph)vhl = (e, + Pq)vy^, (4) 

both sides of the front are at different temperatures and chemical potentials. An impor- 
tant remark here : refers to the quark-gluon plasma side of the sharp front, where we 
allow for an admixture of hadronic matter consisting for simplicity of only pions. Thus 

6, = (1 - /)e; + /e, (5) 

where 6° refers to the pure quark-gluon matter and / is the fraction of hadronic matter 
present in the quark-gluon plasma phase. We will consider / to be a free parameter and 
investigate how various quantities depend on it. It turns out that for the small cut-off 
as well as for the bag model the transition temperature and chemical potential are not 
sensitive to variation of / provided that it is not very large. Fore the large cut-off the 
transition changes completely its nature when / becomes of the order of 20%. The 
thermodynamic quantities (pressure, energy density) refer to the rest frame quantities. 
The index q refers to the quark-gluon plasma side which has an admixture of hadronic 
matter, while the index h refers to the hadronic phase, v refers to the velocity of the gas 
in each phase with respect to the front, 7 is the standard relativistic factor (1/Vl — v'^). 
In addition to the above two equations, baryon number conservation across the front 
leads to 

nhVhlh = riqVq'Jq (6) 

where Uh and Ug refer to the baryon densities on the hadronic and on the quark-gluon 
plasma sides respectively. 

In order to determine the four unknown quantities eq,Pq,Vq and Vh as functions 
of efijPh and Uh we need one more equation. This is given by the requirement of 
non-decreasing entropy which implies in our notation 

ShVhIh > SqVq-fq , (7) 

or, combining (6) and (7) 

. (8) 

rih riq 

As has been shown previously [15], eqs. (3,4,6,7) yield a supercooled quark-gluon 
plasma. The equality sign in (8), corresponding to an adiabatic process, determines 
the minimal amount of supercooling needed to fit the hadronic data without violating 
the second law of thermodynamics. Therefore, we may use the equation 

^ = ^, (9) 
rih riq 
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as the fourth constraint, keeping in mind that the quark-gluon temperature Tq thus 
obtained is maximal. 

Equations (3,4,6) can be solved for the velocities, leading to 




{Ph-Pq){eh + Pq) 
(e/i - eg) (eg + Ph) ' 



(10) 




(11) 



( 
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{eq+Pq){eq+Ph) 

ieh + Ph){eh + Pq)' 



(12) 



Combining now equations (9) and (12) we can calculate numerically the temperature 
Tq and the chemical potential //^. This in turn yields all thermodynamic quantities on 
the quark-gluon side. In the following sections we discuss our main results. 

4. Admixture of Hadronic Matter 

In Fig. 3 we show our results for the case when there is no admixture of hadronic 
matter in the quark-gluon plasma phase. There is a substantial amount of supercooling 
in the quark phase. The transition then takes place to a hadronic gas which has a tem- 
perature T = 200 MeV as suggested by most experiments at CERN using a relativistic 
ion beam. Most results from lattice gauge theories indicate that this temperature is 
higher than the critical temperature of the phase transition. We thus have to conclude 
that the hadronic gas is in a super-heated phase. 

This scenario immediately raises the following question : what happens if other 
mechanisms prevent this strong super-cooling? In particular it has been argued con- 
vincingly by Csernai and Kapusta [21] that super-cooling enhances the formation of 
hadronic bubbles in the quark-gluon plasma phase. We have therefore considered ad- 
mixing hadronic matter to the quark-gluon plasma phase. We expected that the strong 
degree of supercooling would disappear however, we did not expect our results to be so 
sensitive on the value of the cut-off chosen for the equation of state in the quark-gluon 
plasma phase. This is shown in Fig. 4. In this figure we keep the parameters on the 
hadronic side fixed at T = 200 MeV and hb = 300 MeV but we vary the amount of 
hadronic admixture, given by the fraction / as introduced in equation (4), on the quark 
side and we calculate each time the resulting value of the temperature on the other 
side of the sharp front. For simplicity the admixture is in the form of pions only. If 
the cut-off value is large, 5Tc , then there is an abrupt turnover when the admixture 
exceeds 20%. The temperature increase across the front is reduced after this. However 
when the cut-off is taken at 3.4Tc then the picture is completely different : the strong 
supercooling persists anomalously long until the admixture has reached more than 60%, 



4 



in which case it is no longer possible to speak of a "small" admixture. If a bag-model 
equation of state is used the behavior is very similar to the case of the small cut-off. 

In Fig. 5 we show the baryon chemical potential fig on the quark-gluon plasma 
side. The observed behaviour follows the pattern already found for the temperature in 
Fig. 4. As previously, the cut-off in the equation of state is very important. 

5. Velocities 



It is customary to work in the frame where the sharp front separating the two 
phases is at rest. In this frame the quark-gluon plasma moves towards the front with 
a velocity Vq and the hadrons move away from the front with a velocity Vh- If the two 
velocities are different, as is usually the case, one observes a phenomenon of pile up on 
one side, e.g. if the quarks have the largest velocity, the hadrons will pile up with a 
large baryon density at the other side of the front. If Vq < Vh one has a deflagration, if 
Vq > Vh one has a detonation (see [22-24]). 

The velocities can be easily calculated from equations (10) and (11). Our results for 
f g and f ^ as functions of the fraction / of hadronic matter mixed into the quark-gluon 
phase are shown in Figs. 6 and 7 respectively. In Fig. 8 we plot the ratio of v^, to Vq as 
a function of / and in Fig. 9 the relative velocity defined as 

Vh-Vq 

Vrel = z. (13) 

1 - VhVq 

The change in the transition mechanism leads to very sharp variations in Vq and Vh- 
The velocities become unphysical along the curves between the points A and A'. The 
detonation and deflagration branches are clearly separated in Figs. 6 and 7 whereas in 
Figs. 8 and 9 the detonation takes place below the point A and deflagration above the 
point A'. 

As is well known, small perturbations in a medium propagate with the velocity of 
sound in that medium. It is therefore of special interest to compare Vq and Vh with the 
velocity of sound in each phase. The latter is defined as [20] 



2 dP 



(14) 

s/n 



To calculate the speed of sound we write eq. (14) as 

2 _ dP/dT + dP/d^iB ■ d^B/dT 
~ de/dT + de/duB ■ diie/dT ' ^ ' 

where the derivative d^B/dT is taken such that the ratio Sh/uh remains constant. Thus 

d/iB nds/dT - sdn/dT 



dT nds/dfiB — sdn/dfiB 



(16) 
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The velocities of sound in the quark-gluon and in the hadronic phases are shown 
in Figs. 6 and 7 as functions of /. 

The points O, O', A and A' divide the curves in Figs. 6-9 according to the magni- 
tudes of the velocities of sound in the corresponding media. The following inequalities 
hold : 

below O : Vq> Cs,q , Vh < Cs,h 
on OA : Vq < Cs,q , Vh < Cs,h 
on A'O' : Vq > Cs,q , Vh > Cs,h 

above O' :Vq< Cs,q , Vh > Cg^h 

The comparison of Vq and Vh with the velocities of sound in the corresponding 
media shows that the conditions for the stability of the front 

Vq > Cs,q ; Vh < Cs,h (17) 

are fulfilled only along the upper part of the detonation branch. Following the arguments 
presented in Ref. [20] one would conclude that the rest of the curve is unstable and the 
front would disintegrate. Condition (17) is however derived under the assumption of 
a non-decreasing entropy density, i.e. for Sh > Sq. In the present situation we require 
only the entropy per baryon not to decrease which means that Sh < Sq is physically 
acceptable provided the condition expressed in eq. (8) is satisfied. Indeed we find that 
the ratio Sq/sh becomes less than one along the detonation branch (corresponding to 
the part below A' on Figs. 8 and 9) and Sq/ Sh > 1 along the deflagration branch (above 
A'). 

6. Conclusions 

In this paper we have investigated the consequences arising from the existence 
of a sharp front separating the quark-gluon plasma phase from the hadronic phase. 
The equations of continuity for energy, momentum and baryon number as well as the 
constraint following from the second law of thermodynamics (entropy increase) , require 
a substantial amount of supercooling to occur in the quark-gluon plasma before the 
transition can take place. This transition then leads to a highly super-heated hadronic 
gas. We have also investigated what happens if the quark- gluon plasma is not pure 
but contains an admixture of hadronic matter. It is known that supercooling is limited 
when "impurities" are present in the super-cooled phase. We found that from the 
moment when the admixture of hadronic matter is more than about 20% of the total 
(the fraction refers to the energy content, not to the particle content) then the transition 
does not take place via a detonation type of wave but via a deflagration one. In this 
case the transition is smoother and the temperature jump across the sharp front is not 
as pronounced as in the absence of admixture of hadronic matter. 

The crossover region between the two solution shows a remarkable S-shape struc- 
ture. The multiple valued Tq and in that region indicate that in the given range of 
/ a chaotic turbulent type of transition could take place. 
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In summary, we have considered a model based on the existence of a sharp front 
separating the hadronic from the quark-gluon phase. We have arrived at the conclusion 
that in this situation the transition can only proceed if a substantial amount of super- 
cooling takes place in the quark-gluon plasma phase if the transition is between pure 
phases. An admixture of hadronic material in the quark-gluon plasma phase makes the 
transition proceed in a much smoother way. 
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Figure Captions: 

1: Energy-density as a function of temperature at fixed hb ~ 0.3 0.3 GeV. The long- 
dashed fine corresponds to the equation of state of the M.I.T.-bag with bag constant 
B=0.2 GeV/fm^, the dashed fine is the cut-off model with a value 3.4Tc for the cut- 
off while for the long dash-dotted line the cut-off is 5Tc. The solid line corresponds 
to the energy-density as calculated from the hadronic gas. The critical temperature 
is indicated by the dotted line at Tc = 150 MeV. 

2: The phase diagram calculated in the cut-off model and in the M.I.T. bag model. 
At fiB = the bag constant is chosen such as to reproduce the same critical 
temperature as the cut-off model. The value of the cut-off is taken to be 5Tc 
independent of fiB- 

3: The energy-density of the hadronic gas (lower curve) and of the quark-gluon plasma 
(upper curve) as a function of temperature. Super-cooling and super-heating are 
indicated by the dashed line. The transition from the super-cooled quark-gluon 
plasma phase to the hadronic gas phase is indicated by the long-dashed line. 

4: The quark temperature Tg versus the fraction of hadronic admixture in the quark- 
gluon plasma phase. The solid line corresponds to a cut-off value of 5Tc, the dashed 
line to 3.4Tc. The dash-dotted line is for the equation of state of the M.I.T.-bag 

model. 

5: The baryon chemical potential in the quark-gluon plasma phase versus the fraction 
of hadronic admixture in the quark-gluon plasma phase. The solid line corresponds 
to a cut-off value of 5Tc, the dashed line to 3.4Tc. The dash-dotted line is for the 
equation of state of the M.I.T.-bag model. 

6: The velocity squared of the quark matter (full line). The value of the cut-off pa- 
rameter in the equation of state is 5Tc, / is the admixture fraction. The dashed line 
is the velocity of sound squared in quark matter. See text for further explanations. 

7: The velocity squared of the hadronic gas (full line). The value of the cut-off param- 
eter in the equation of state is 5Tc, / is the admixture fraction. The dashed line is 
the velocity of sound in the hadronic matter. See text for further explanations. 

8: Ratio of the hadronic gas velocity over the quark phase velocity as a function of the 
admixture fraction f. The value of the cut-off parameter in the equation of state is 
5Tc. If the ratio is less than one, the transition proceeds via a detonation wave, if 
it is above one, it proceeds via a deflagration wave. Between the points A and A' 
the velocities are unphysical. See text for further explanations. 

9: The difference between the two velocities. The value of the cut-off parameter in 
the equation of state is 5Tc, / is the admixture fraction. Between the points A and 
A' the velocities are unphysical See text for further explanations. 
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